Dissipation-induced pure Gaussian state 
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This paper provides some necessary and sufficient conditions for a general Markovian Gaussian 
master equation to have a unique pure steady state. The conditions are described by simple matrix 
equations, thus the so-called environment engineering problem for pure Gaussian state preparation 
can be straightforwardly dealt with in the linear algebraic framework. In fact, based on one of those 
conditions, for an arbitrary given pure Gaussian state, we obtain a complete parameterization of 
the Gaussian master equation having that state as a unique steady state; this leads to a systematic 
procedure for engineering a desired dissipative system. We demonstrate some examples including 
Gaussian cluster states. 

PACS numbers: 03.65.Yz, 42.50.-p, 42.50.Dv 



I. INTRODUCTION 

Preparing a desired pure state, particularly under in- 
fluence of dissipation, is clearly a most important sub- 
ject in quantum information technologies. To tackle this 
problem, other than some well-acknowledged strategies 
such as quantum error correction, recently we have a 
totally different method that rather utilizes dissipation 
0, S II 1, II S S 0, i M, HO, M O ■ The basic idea is 
to engineer a dissipative system so that the system state 
p(t) governed by the corresponding Markovian master 
equation 
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must evolve towards a desired pure state: p(t) — > |</>)(</>| 
as t — > oo. Here, H is the system Hamiltonian and L k 
(k = 1, . . . , m) is the dissipative channel that represents 
the coupling between the system and the fc-th environ- 
ment mode. The main advantage of this environment en- 
gineering approach is that the dissipation-induced state 
\4>) is robust against any perturbation and thus may serve 
as a desired state, e.g., an entangled state for quantum 
computation Q and quantum repeater [l2j]- Therefore, a 
complete characterization of the master equation having 
a unique pure steady state should be of great use, and ac- 
tually in the finite-dimensional case it was given by Kraus 
et al. [4]. In particular, they showed that some useful 
pure states including cluster states can be prepared by 
quasi-local dissipative process, i.e., dissipative channels 
that act only on a small number of subsystems. 

In this research direction the infinite-dimensional coun- 
terpart of the above-mentioned state preparation method 
should be explored. In particular, Gaussian states con- 
stitute a wide and important class of quantum states, 
which serve as the basis for various continuous- variable 



(CV) quantum information processing [Tj, The con- 
tribution of this paper is to provide some necessary and 
sufficient conditions for the master equation ([1]) to have a 
unique pure steady state, when H and L k are of general 
form for p{t) to be Gaussian for all t. The conditions 
are described by simple matrix equations, thus they can 
be properly applied to the above-mentioned environment 
engineering problem for pure Gaussian state preparation. 
Actually, one of those conditions enables us to obtain a 
complete parameterization of the Gaussian master equa- 
tion that uniquely has a pure steady state. This leads 
to a systematic procedure for constructing a dissipative 
system deterministically yielding a desired pure Gaus- 
sian state. We provide some examples of dissipation- 
induced states including the so-called Gaussian cluster 
states [lH, [I?], EH , which are known as essential resources 
for the CV one-way quantum computing [l9| , with focus- 
ing on how they can actually be prepared by quasi-local 
dissipative process. 



II. GAUSSIAN DISSIPATIVE SYSTEMS 

We here provide the phase space representation of the 
general Gaussian dissipative system with n-degrees of 
freedom, which is subjected to the master equation ([1}. 
Let (<ji, pi) be the canonical conjugate pair of the i-ih 
subsystem. It then follows from the canonical commu- 
tation relation [g,-, pj] = iSij that the vector of system 
variables x := (<?i, . . . , q n ,Pi, ■ ■ ■ ,Pn) T satisfies 
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where /„ denotes the nxn identity matrix and T the ma- 
trix transpose. A Gaussian state is completely character- 
ized by only the mean vector (x) and the (symmetrized) 
covariance matrix V = (AxAx T + (AxAx T ) T )/2, Ax = 
x — (x). Here the mean (X) = Tr (A/5), with p the cor- 
responding Gaussian density operator, is taken element- 
wise; e.g., (x) = ((<?i), . . . , (p n )) T ■ The uncertainty rela- 
tion is represented by the matrix inequality V+iT,/2 > 0, 
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thus V > 20] . Note also that the purity of a Gaussian 
state is simply given by P = Tr (p 2 ) = l/^/2 2n det(F). 
Now consider Eq. ^ with H and given by 

H = ^x T Gx, L k = cjx, 

where G = G T G R 2 " x2 " and c fe G C 2 " are the param- 
eter matrix and vector specifying the dissipative system. 
Then, the time-evolution of (x) and V read 

"-f f -^ + VAT + Dl (2, 

where A = E[G + 9f(C*C)] and £> = E9?(GtG)E T with 
G = (ci,...,c m ) T G C mx2n (3? and 3 denote the real 
and imaginary parts, respectively). If p(0) is Gaussian, 
then p(t) is always Gaussian with mean (x(t)) and co- 
variance matrix V(t). See [2l[ for detailed description. 

Now suppose that A is Hurwitz; i.e., all the eigenvalues 
of A have negative real parts. This is equivalent to that 
the system has a unique steady state; it is Gaussian with 
mean (x(oo)) = and covariance matrix V s that is a 
unique solution to the following matrix equation: 

AV S + V S A T + D = 0. (3) 

Note that V s can be explicitly represented as 

/>oo 

V s = / e At De ATt dt. (4) 
Jo 

The purpose of this paper is, as mentioned before, to 
fully characterize a Gaussian master equation that has 
a unique pure steady state, and this is now expressed in 
terms of Eq. (H]) by 2 2n det(14) = 1. However, clearly 
this is not useful. In the next section we give a much 
simpler and explicit version of such a characterization. 

III. THE DISSIPATION-INDUCED PURE 
GAUSSIAN STATES 

A. Pure steady state condition 

Here we address our first main result: 
Theorem 1: Suppose that Eq. ^ has a unique solution 
V s . Then the following three conditions are equivalent: 

(i) The system has a unique pure steady state with 
covariance matrix V^. 

(ii) V s satisfies the following matrix equations: 

(V s + ^e)g t =0, (5) 
EGV a + V s GYi T = 0. (6) 

(iii) The following matrix equation holds: 

KY,C T = 0, (7) 
where K G C 2nmx2 ™ is defined by 
K = (G T , GE T G T , • ■ ■ , (GS T ) 2 "- 1 C T ) T . (8) 



Furthermore, when the above equivalent conditions are 
satisfied, V s is represented by 

V s = is T S(^tx)[SR(Xtif)]- 1 . (9) 

To prove this theorem we need the following lemma. 
Lemma 1: Suppose that Eq. ([3]) has a unique solution. 
Then rank(A' T ) = 2n, where 

K = (C T , GS T C T , . . . , (GE T ) 2 "- 1 C' T ) T G R 4 ™ mx2 ™ 

with C = (3?(G) T ,3(G) T ) T G R 2mx2 «. 

Proof of Lemma 1: From the assumption, Eq. ([3]) 
has a unique solution (j4j. Suppose there exists a vector 
£ G R 2 ™ such that = 0. Then, noting that A and D 
are written by A = E[G + G T EG] and D = EG T GE T , 
we have £ T E T y s E£ = / °° \\CH T e^^fdt = 0, and 
this is contradiction to V s > 0. Thus rank(X T ) = 2n. ■ 

Proof of Theorem 1: The proof is divided into four 
steps. 

1. (i) => (ii). From the Williamson theorem [22| a co- 
variance matrix V corresponding to a pure state is ex- 
pressed by V = SS T 1 2 with S a symplectic matrix [23| . 
Thus substituting V s = SS T /2 for Eq. ©, multiplying 
it by S^ 1 from the left and by S~ T from the right, and 
finally using this equation twice to erase G, we have 

(g; t -eg 1 ' t )(g; t -egi t ) t 

+ (CJ + eg; t )(g 2 ,t + eg; t ) t = o, 

where C' r = U(C)S and C[ = <S(C)S. Note SY,S T = E. 
Thus G, ,T - EGj' T = G; T + EG; T = 0, which immedi- 
ately leads to V s $t(C) T = E3(G) T /2 and V S ${C) T = 
— E3?(G) T /2. From these equations we obtain Eq. (|5|). 
Combining this with Eq. (|3|) yields Eq. ([6]) . 

2. (ii) ^> (iii). Multiplying Eq. flg) by EG from the 
left and using Eq. ©, we get (V s + iE/2)GE T G T = 0. 
Repeating this manipulation yields 

(V s + i^/2)K T = 0, (10) 

with K defined in Eq. ©. This readily leads to K(V S ± 
iY,/2)K T = 0, thus KY,K T = 0. The Cayley-Hamilton 
theorem implies Eq. (|7]l. 

3. Derivation of Eq. ([9]). Define the 4nm x 2n matrix 
K' = ($t(K T ), 3(X T )) T , then the real and imaginary 
parts of Eq. (JTUJ) are summarized in a single equation as 

YV S K' T = K' T E/2. 

Now, from the assumption, we can use Lemma 1 and 
find rank(iT' T ) = rank(X T ) = 2n. Then multiplying the 
above equation by K' (K' T K')^ 1 (the generalized inverse 
matrix of K' T ) from the right, we have 

EV S = (KjKi - KjK r )(KjK r + Kj K^ 1 /2, 

where K r = $t(K) and K { = %(K). This is just Eq. ©. 



3 



Now let us verify that Eq. ([9]) is symmetric. Noting 
that Eq. ([7]) is equivalent to KEK T = 0, we have 



K 

K* 



E{K 



-AT 



K 

-K* 



E T (K T , Jift). (11) 



Then multiplying this equation by {K\ K T ) from the 
left and by (A^, K T ) T from the right, we obtain 

K(A t A)SS(A t A) = -^(K^K)E T ^(K^K), 

thus V s = Vj. 

4. (iii) =>■ (i). First, to show that the state is pure, 
we use the fact [25j that, for a pure Gaussian state, the 
corresponding covariance matrix V satisfies 



EVEV 



-1/4. 



(12) 



Now multiply Eq. (JTTJ) by {K\ K t ) from the left and 
by (At, — K T ) T from the right, then we have 

^.{K^ A)S3J(A t A) = Q?( J ftr t X)I] T Q?( J ft: t X). 

This equation readily implies that V s given by Eq. © 
satisfies Eq. (IT21) . hence the corresponding state is pure. 

Next, let us show that Eq. @ satisfies Eq. ©. Note 
that Eq. © is rewritten as V s $t{K^K) = -EQ(K^K)/2, 
Also from Eq. (I2J) we have ^(A^A) = E5R(Atii')/2. 
These two equations yield 

(F s - iS/2)A + A = 0, 

which is equivalent to (F s - i£/2)At = 0, thus Eq. ([TO]). 
This implies that Eq. © holds. Moreover, Eq. ([TP)) 
leads to (V S GE T + i£G£ T / 2 )A T = and (EGV a + 
iEGE/2)K = 0; from these equations we have (EGV S + 
V S GE T )K T = 0, and as now EGV a + V S GE T is real and 
rank(A'~) = 2n, we obtain Eq. ©. As a result, V s satis- 
fies Eqs. ([5]) and (|6|), but these two equations correspond 
to a specific decomposition of Eq. That is, V s is the 
solution to Eq. ©. ■ 
We give an interpretation of Theorem 1 . Let p s be the 
pure density operator corresponding to the covariance 
matrix V s . Then Eqs. ([5]) and (JBJ) are equivalent to 

2L k p s L{-L{L k p s -p s L{L k = Vfc, [H, p s ] = 0, (13) 

respectively. The former condition is further equiva- 
lent to that Lk\4> s ) is parallel to \4> a ) for all k, where 
p s = \(p s )(<p s \. This means that p s is the so-called dark 
state; that is, Eqs. (J5J) and (JB]) are the phase space repre- 
sentation of the condition for the state to be dark. More- 
over, the uniqueness of p s allows us to erase itself in Eq. 
(fTB")) and obtain a single equation with respect to H and 
Lfc. The phase space representation of this equation is 
no more than Eq. ([?])• It should be maintained that p s is 
explicitly represented in a directly computable form (|9|). 

Next let us discuss how to use Theorem 1, particularly 
for the purpose of environment engineering. First, note 
that Eq. (J7J) depends only on the system matrices G and 



C; thus the condition (iii) should be applied, for a given 
specific system configuration, to find the system param- 
eters such that the corresponding master equation has a 
unique pure steady state. On the other hand, the con- 
dition (ii) explicitly contains V^; this means that we can 
characterize the structure of a dissipative system such 
that a desired pure state with covariance matrix V s is 
generated by that dissipative process. Later we provide 
a modification of the condition (ii) that can be more suit- 
ably used to find such a dissipative system. 



B. Examples 

We here give two examples to explain how the theorem 
is used. 

Example 1: Single OPO. Let us first study an ideal 
optical parametric oscillator (OPO), which couples with 
a vacuum field through one of the end-mirrors. The 
Hamiltonian is described in terms of the annihilation 
operator a = (q + ip) /y/2 and the creation operator 

= (q — ip)/\/2 as H = i(ea) 2 — e*a 2 )/4, where e de- 
notes effective complex pump intensity proportional to 
X*- 2 - 1 coefficient of the nonlinear crystal. Also the cou- 
pling operator is given by L = ^fna with k the damping 
rate of the cavity. The corresponding matrices (G, G) 
then take the following form: 

G -{ 3(e) -»(e) J ' 6 " V 2 (1 ' l) - 

From Theorem 1, the steady state of this system becomes 
pure if and only if 



KEG" = 



GSG T 
GSGEG 1 







= 0. 



That is, the dissipative process brought about by the 
coupling to the outer vacuum field must introduce deco- 
herence to the intra-cavity state as long as it is squeezed 
(e 7^ 0). Actually, when e = the steady covariance 
matrix ^ turns out to be I2/2, In conclusion, a single- 
mode intra-cavity state can become pure only when it is 
the trivial vacuum (or coherent) state. I 
Example 2: Cascaded OPOs. We next consider the 
two- mode OPOs shown in Fig. 1 (a), where the OPOs 
are connected through a unidirectional optical field. This 
kind of cascaded system plays an important role in build- 
ing a quantum information network; e.g., entanglement 
distribution was discussed in Q. The physical setup of 
each OPO is the same as before, i.e., Hj — iej(aj 2 — aj)/4 

and Lj = ^/naj (j = 1,2) with d\ and 02 each cavity 
modes. For simplicity we here set the squeezing effective- 
ness of each cavity to be real; ej £ R. From the theory of 
cascaded systems [2^, [27j , the Hamiltonian and the cou- 
pling operator of the whole two-mode Gaussian system 
are respectively given by 

H = H 1 +H 2 + (L\U - L\L 2 )/2i, L = L X + L 2 , 
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hence the corresponding system matrices read 



OPOl 



G 





( 





ei 


— K 










K 


£2 




I £1 


K 












£2 









c 



(1, 1, i, i). 



Then A = E(G + %(C^C)) has eigenvalues -(« ± e,-)/2, 
thus it is Hurwitz when |ej| < k (j — 1, 2). Equivalently, 
when both the OPOs are below threshold, the steady 
state is unique. Then from the condition (iii) of Theo- 
rem 1 the steady state is pure if and only if 



KY.C 1 = 



CEC T 
GEG£G T 
C(£G) 2 £C T 
V G(SG) 3 EG T J 



i«(ei + e 2 ) 



= 0, 



where A = e\ + e| — 6162 — 3k 2 7^ 0. Therefore, the system 
should be engineered to satisfy £1 + 62 = to have a 
unique pure steady state. From Eq. the corresponding 
covariance matrix is given by 
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FIG. 1: (Color online) (a) Schematic diagram of the cascaded 
OPOs. (b) Time evolutions of the fidelity F(t) and the purity 
Pit) of the bipartite cavity state. The parameters are k = 6.0 
MHz and e — 4.8 MHz. We take the initial Gaussian states 
with mean zero and covariance matrices V(0) = 1/2 (top 
line), V(0) = I (middle), and V(0) = 21 (bottom). Note 
due to (x(t)) = Vt, the fidelity is now of the form F(t) = 
1 / y/det(V(t) + V s )- The entanglement achieved in this case 
is E(V S ) = 1.0986. 



Uniqueness condition for a given Gaussian pure 
steady state 



v. = - 
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K/g+ 


£/3+ 








t/g+ 


K /g+ , 



(14) 



where g± :— k ± e and e := ei = —62. Thus the steady 
state is a nontrivial entangled pure state other than the 
vacuum when < \e\ < k. Actually the logarithmic 
negativity [28] , which is a convenient computable measure 
for entanglement, takes a positive value; 



We again stress that this dissipation-induced entangled 
state is guaranteed to be highly robust. That is, any ini- 
tial state /3(0) converges into that entangled state. Such 
robustness can be clearly observed from Fig. 1 (b) that 
demonstrates the time-evolutions of the fidelity F(t) = 
Tr (p{t)p s ) and the purity P{t) = Tr (p(t) 2 ) with several 
initial states. ■ 



IV. ENVIRONMENT ENGINEERING FOR 
PURE GAUSSIAN STATE PREPARATION 



In this section we first address a modification of the 
condition (ii) of Theorem 1 with different assumption, 
which is more suited to the concept of environment engi- 
neering. In fact, the result allows us to obtain a general 
procedure for synthesizing a dissipative system whose 
steady state is uniquely a desired pure Gaussian state. 
We close this section with some examples. 



Theorem 2: Let V s be a covariance matrix correspond- 
ing to a pure Gaussian state. Then, this is a unique 
steady state of the system if and only if 



ker (v s + ^Y^j — range(K 1 



(15) 



where K is defined in Eq. (J5J). 

Proof: We begin with the necessary part. In general, 
for a n-mode Gaussian pure state, V s + iS/2 has a n- 
dimensional kernel (23|, while now from the assumption 
Eq. (HD} holds, implying range(i^ T ) C kei(V s + iS/2). 
Hence, showing rank(i^ T ) = n completes the proof. 
Now, as shown in the proof of Theorem 1, from the as- 
sumption we have KY,K T = 0, implying rank(.ftT T ) < n. 
On the other hand, the uniqueness of the steady state 
allows us to apply Lemma 1 to get rank(^ T ) = 2n, thus 
rank(if T ) > n. As a result we obtain rank(i\T T ) = n. 

Let us next move to the sufficiency part. First we 
show that V s satisfying Eq. (|T5|) is a solution of Eqs. 
© and ©, thus that of Eq. ©. Note that now Eq. © 
apparently holds. As seen from the last part of the proof 
of Theorem 1, Eq. $T5§ implies (ZGV S + V S GZ T )K T = 0, 
thus we need rank(i^' T ) = 2n to derive Eq. ([6]). To show 
this, assume there exists £ e R 2 ™ such that K' £ = 0. 
This leads to K£ = 0. Now, from Eq. ([T5) with V s 
pure, we have KHK T — and rank(X T ) = n, thus 
ker (if) = range(Eif T ). Then, as £ € range(Si4T T ), we 
can write £ = ZK 1 a, 3a e C 2nm . This leads to = 
U[(V„ +iZ/2)K T a] = V S E T £, thus contradiction to V s > 
0. Consequently, we have rank(if' T ) = In. 

Second, to complete the sufficiency part, we show that 
V s is unique, which is equivalent to that A is Hurwitz. 
Suppose A T ^ = A£ with £ e C 2 " and A e C. Multiplying 
Eq. ^ by ^ from the left and by £ from the right, we 
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obtain »(A) = -^Df/W^f)- Note in general D > 
and V, > 0. But now > 0, because = leads 

to X'I] T £ = and this is contradiction to rank(if' T ) = 
2n. Thus 3f (A) < and this means that A is Hurwitz. ■ 
The point of this result is that, while in Theorem 1 
the steady state is assumed to be unique, we here as- 
sume only that a given state is pure, without assuming 
its uniqueness. Nonetheless, that state is guaranteed to 
be a unique steady state if the condition (fPoj) is satisfied. 
Note again, as seen from the last part of the above proof, 
that the uniqueness is ensured by rank(if T ) = n, which 
leads to the Hurwitz property of the matrix A. That is, 
we do not need to check whether A is Hurwitz. 



B. Complete parameterization of the dissipative 
system 

Based on the result shown above, we here provide a 
complete parameterization of the Gaussian dissipative 
system that uniquely has a pure steady state. This then 
leads to an explicit procedure for engineering a desired 
Gaussian dissipative system. 

We begin with the fact that any covariance matrix V s 
corresponding to a pure Gaussian state has the following 
general representation [l8|, [24[ : 

Vs = -SS T , S = I ^ y _ 1/2 Y \ /2 J , (16) 

where X and Y are nxn real symmetric and real positive 
definite matrices (i.e., Y = Y T > 0), respectively. Note 
that S is symplectic. With this representation we have 

where we defined Z = X + iY . It was shown in [l8| that 
the symmetric matrix Z is useful in graphical calculus for 
several Gaussian pure states. Because (I, Z) is clearly of 
rank n, we have 

ker (v s + ~X!^ = range ( , . 



Again, P (complex), R (real symmetric), and T (real 
skew) are the parameter matrices. We now have a rea- 
sonable procedure for environment engineering for pure 
Gaussian state preparation; that is, the procedure is sim- 
ply to choose the matrices P, R, and T so that both 
the dissipative channels — cjx (k = l,...,m) and 



Hence, Eq. (|T5j) is equivalent to 

range ( ~/ ) = range(G T , G£ T C T , . . . , (GE T ) 2 "~ 1 C T ). 

To satisfy this condition, it is necessary that range(G T ) 
is included in range(— Z, I) T and that range(— Z, I) T is 
invariant under GS T . These conditions are respectively 
represented by 

C t =(~/)p, GE T ("/) = ("/)q, (17) 

where P and Q are n x m and nxn complex matrices. 
From the above equations, K T is represented by 

K T = (~f^j (P, QP,..., Q^-'P). 

Consequently, the necessary and sufficient condition for 
range(if T ) to be identical to range(— Z, T) T is that there 
exist P and Q satisfying Eq. (fT7| and the rank condition 

rank(P, QP, . . . , Q n ~ l P) = n. (18) 

Now let us write G in the 2x2 block matrix form 
G = (Gi, G 2 ; Gj, G3), where G\ and G3 are real nxn 
symmetric matrices and G 2 a real nxn matrix. Then 
the latter equation in Eq. (fTTj) leads to 

Gi + G 2 X + XG 2 + XG 3 X - YG 3 Y = 0, 
(G 2 + XG 3 )Y + Y(G 2 + XG 3 ) T = 0. 

The second equation is equivalent to that there exists a 
real skew symmetric matrix T (i.e., T+r T = 0) satisfying 
(G 2 + XG 3 )Y = T. Hence, by writing R = G 3 , we find 
that G2 is expressed as G2 = —XR + TV -1 . In this 
representation, Q is of the form Q = —iRY — Y~ x r T . 
To conclude, we obtain the complete parameterization of 
G and G as follows: 



(19) 

I 

the Hamiltonian H — x T Gx/2 with the system ma- 
trices given in Eq. (|19[) have desired structures such 
as quasi-locality, while, at the same time, P and Q = 
-iRY - F~ 1 r T satisfy the rank condition (fI5|). 

Here we give some remarks. 

(i) There always exists the pair of (P, R, T) satisfying 



C = P T (-Z, I), G 



XRX + YRY -TY^X - XY- 1 ^ -XR + TY 
-RX + Y- 1 ^ R 
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the rank condition Eq. (jl8|) . if no restriction is imposed 
on those matrices; this means that, for any pure Gaussian 
state, there always exists a Gaussian dissipative system 
for which that state is the unique steady state. 

(ii) The most simple system may be such that P is of 
rank n. In this case we can set R = T = 0, implying that 
the system does not need a nontrivial Hamiltonian but 
drives the state only by dissipation. This kind of system 
is called the purely dissipative system. However, we are 
often in the situation where only m < n dissipative chan- 
nels can be implemented in reality, due to some reasons 
related to physical constraints. In this case, the matrix 
Q needs to be of rank at least n — m, meaning that we 
must add a nontrivial Hamiltonian. 

(iii) As stated in quasi-locality is indeed essential 
since otherwise it would be experimentally hard to realize 
such a dissipative system. It should be noticed that, be- 
cause H is quadratic, it can be always decomposed into 
the sum of quasi-local Hamiltonians acting on at most 
two nodes, although those interactions between the nodes 
do not necessarily have the structure of a target entan- 
gled state; in fact, it will be shown in Example 5 that, to 
generate a chain-type cluster state, a Hamiltonian having 
a ring-type interaction is added. That is, while in Gaus- 
sian case the quasi-locality issue appears only in the part 
of dissipative channel, this does not mean that the com- 
plementary Hamiltonian can readily be implemented. 



C. Examples 

Example 3: General CV cluster and %-graph states. 
Menicucci et al. developed in [l8[ a unified graphical cal- 
culus for all pure Gaussian states in terms of the matrix 
Z = X + iY. One of the important results is that the 
so-called canonical CV cluster state, which can be gener- 
ated by first squeezing the momentum quadrature of all 
modes and then applying the controlled Z operations to 
the modes according to the graph of the cluster, can be 
generally represented by 



Z 



X + ie~ 2r L 



(20) 



where X corresponds to the symmetric adjacency matrix 
representing the graph structure of the cluster state and 
r the squeezing parameter. For this state, for instance 
setting P = I and Q = in Eq. fl"?)) gives a desired 
purely dissipative system with channels 



(-x 



-2r 



I,I)X, 



which have the same structure as that of the target clus- 
ter state. Therefore, if each node of the graph is con- 
nected to at most three adjacency nodes, each dissipa- 
tive channel acts on at most three modes too, i.e., it is 
quasi-local. 



Another important state discussed in [18| is the %- 
graph state; this state is generated by applying the uni- 
tary transformation U = exp(— iHt/K) with Hamiltonian 



H = ihn 



3,k 



W jk {a)a\ 



ajd k ) 



(21) 



to the vacuum states |0) , where W = (Wjk) is the 
real symmetric matrix representing the graph. Note that 
H is the sum of the two- mode squeezing Hamiltonians. 
The corresponding covariance matrix is then given by Eq. 
dUD with X = and Y = e~ 2aW , hence 



Z 



-2aW 



(22) 



where a = 2nt. Unlike the CV cluster state representa- 
tion (f2"0"l) , Z does not necessarily reflect the graph struc- 
ture of the state. However, for instance when W is 
self-inverse, i.e., W 2 = /, we have Z = icosh(2a)7 — 
ismh(2a)W. Thus, in this case choosing P = I gives a 
desired purely dissipative system acting on the nodes in 
the same manner as the Hamiltonian (|2ip . B 
Example 4- Two-mode squeezed state. In the Gaus- 
sian formulation we are often interested in the two-mode 
squeezed state, as it approximates the so-called EPR 
state. This is represented as a H-graph state with the 
Hamiltonian (|2"Tj) given by 



W = 



1 

1 



H = 2ihn(a\a\ — aid^) 



Then the system matrices are X = and 
Y 



-2aW 



cosh(2a) — sinh(2a) 
— sinh(2a) cosh(2a) 



The corresponding covariance matrix (|16[) is given by 



V '=2 



V 



cosh(2a) 


sinh(2a) 








sinh(2a) 


cosh(2a) 














cosh(2a) 


— sinh(2a) 








— sinh(2a) 


cosh(2a) , 



Let us begin with constructing a purely dissipative sys- 
tem whose unique steady state is the above two-mode 
squeezed state; in this case, we only need to specify a 
2x2 full-rank matrix P to satisfy the rank condition 
(|18p . In particular, let us here choose 



P = 



i cosh(a) i sinh(a) 
isinh(a) zcosh(a) 



(23) 



which is clearly of full rank. Then, C — P T (~Z,I) 
(ci,c 2 ) T is given by 



ci 



(p, v, in, -w) T , 



C2 



(v, fi, -iv, ifi) T , (24) 



where /i = cosh(a) and v = — sinh(a). As a result, the 
master equation describing this purely dissipative process 
is given by 



(25) 



fc=l,2 
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where L\ = /iai + va\ and Li — [iai + va\. A possible 
physical realization of this system in a pair of atomic 
ensembles was discussed in pjj, LLJJ • 

Next let us discuss the case where we are allowed to 
implement only one dissipative channel. As an example 
we take P — (i cosh(a), i sinh(a)) T , which corresponds to 



the first column vector of Eq. ([23)) . This yields C = c[ 
in Eq. (|24[) . thus the dissipative channel is L\ = fiai + 



va\. We now need to specify a valid Q to satisfy the 
rank condition (|18[) ; in particular let us take a specific 
Hamiltonian matrix G in Eq. (fT9|) with R — diag{0, 1} 
and r = 0, then this leads to 



q = —iRY - y _1 r T = 



o o 

i sinh(2a) —i cosh(2a) 



It is easily verified that (P, QP) is of full rank, hence 
the requirement is satisfied. In this case the Hamiltonian 
H = x T Gx/2 is given by 



H = (<7i sinh(2a) — q 2 cosh(2a)) 2 



This Hamiltonian and the dissipative channel L = fiai + 
va\ construct the desired dissipative system. ■ 
Example 5: 1- dimensional harmonic chain. As a typ- 
ical cluster state let us take a 1-dimensional (equally 
weighted) harmonic chain, particularly in the case of 
four-mode cluster just for simplicity. Within the for- 
malism of the canonical CV cluster state generation, the 
adjacency matrix X and the graph matrix (|20j) are re- 
spectively given by 



X 



A desired purely dissipative system is readily obtained 
by taking P = J 4 in Eq. (17]), i.e., C = (-Z, 7); this 
means that the four dissipative channels are given by 





1 








( : 





1 




H 






1 











1 


J 






U 


— (— ie 


' 2r qi 4 


■Pi)- 


82, 






U 


= 


- (-ie 




HP2)- 


83, 




u 


= -82 - 


- (-ie 


-2r~ 
<73 ' 




84, 




u 


= -83 - 


- (-ie' 


-2rz 
94 - 






(26) 



Each channel acts on at most three nodes, thus they 
are quasi-local; Fig. 2 (a) depicts the structure of this 
environment-system interaction. Note from the above 
discussion that the general 1-dimensional harmonic chain 
can also be generated by purely dissipative process with 
their channels acting on at most three nodes. The finite 
dimensional counterpart to this result is found in [4]. 

Now we have a natural question; can the chain state be 
generated by a quasi-local dissipative process acting on at 
most two adjacency nodes? If this is true, this means that 
engineering the dissipative environment becomes easier, 
apart from that we clearly need an additional Hamilto- 
nian. Again let us consider the case of four-mode chain 



and take P — (1, 0, 0, 0) T , implying that the system has 
one dissipative channel L\ in Eq. (|26l) . Note this acts on 
only two nodes. To determine the Hamiltonian we have 
some freedom, but let us take in Eq. (TT5)) 



R = X- 1 = 



10-1 
10 
1 
-10 10 



r = o. 



This leads to Q = -iRY - y- 1 r T = -ie^X' 1 , and it 
is readily verified that (P, QP, Q 2 P, Q 3 P) is of full rank. 
Hence, we obtain the positive answer to the question 
raised above. The matrix G is now of the form 



G 



X 



-2r Y-l _J 

x- 1 



e - r X 
-I 



The first (1,1) block matrix indicates that the Hamilto- 
nian H = x T Gx/2 has a ring-type structure where the 1- 
2, 2-3, 3-4, and 4-1 nodes are connected with each other; 
see the remark (iii) in the previous subsection. Therefore 
it is concluded that the four-mode harmonic chain state 
can be generated by the dissipative process where both 
the dissipative channel and the Hamiltonian quasi-locally 
act on only two adjacency nodes; the structure of these 
interactions is shown in Fig. 2 (b). Implementation of 
this system should be easier than that of the purely dis- 
sipative system obtained above, which acts on the nodes 
through the channels (|26|) . ■ 



• 3 




FIG. 2: (Color online) (a) Dissipative channels for generating 
the harmonic chain state, (b) A combination of the dissipative 
channel (arrow) and the ring-type Hamiltonian (dotted circle) 
can generate the harmonic chain state. 



V. CONCLUSION 

In this paper, we have derived the general necessary 
and sufficient conditions for a Gaussian dissipative sys- 
tem to have a unique pure steady state. In particular, we 
have provided a complete parameterization of the Gaus- 
sian dissipative system satisfying the requirement; this 
leads to an explicit procedure for engineering a Gaus- 
sian dissipative system whose steady state is uniquely a 
desired pure state. 

An important open question is how to construct, in 
a general framework, a practical dissipative system sat- 
isfying a quasi-locality constraint. Although Example 5 
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demonstrated the feasibility in engineering such a quasi- 
local dissipative system, that construction method is a 
heuristic one. Note that in the finite-dimensional case 
the same problem remains open 

We can take two approaches to dealing with the prob- 
lem. The first one is suggested by the recent work by 
Ticozzi and Viola 

.13] that, in the finite-dimensional case, 
characterizes a pure state generated by a dissipative sys- 
tem having a fixed quasi- locality constraint; it was then 
shown that, based on this result, by switching dissipa- 
tive channels through output feedback control [2l|, we 
are enabled to construct a desired quasi-local dissipative 
system. It is expected that this idea works in our case as 
well. Actually, it was shown in Q that, with the use of a 
similar switching method, some Gaussian cluster states 
can be generated in four-mode atomic ensembles trapped 
in a ring cavity. 

The other approach uses the schematic of quantum 
state transfer. The basic idea is as follows. First, an n- 
mode target (entangled) Gaussian state of light fields is 



produced, and then, those fields are independently cou- 
pled to n identical bosonic systems in a dissipative way; 
then, as a result, that state is deterministically trans- 
ferred to the systems. In this scheme, we need the strong 
assumption that all the nodes are accessible and indepen- 
dently couple to the fields, but the interactions between 
the systems and the environment channels are all local. 
Hence, from the fact that a number of pure Gaussian clus- 
ter state of light fields can be produced efficiently with 
the use of some beam splitters and OPOs [MITEEl, this 
state-transfer-based approach is expected to be more rea- 
sonable. In fact, the dissipation-induced states shown in 
[H G3 EH are generated using this technique. Further- 
more, one of the authors recently has developed a general 
theory of this scheme in terms of a quantum stochastic 
differential equation [2PJ. 

It would be worth to further examine the above two 
approaches to make them more useful for engineering a 
practical dissipative system. 
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